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Abstract 

The present note contains the text of lectures discussing the problem of universal- 
ity in fully developed turbulence. After a brief description of Kolmogorov's 1941 
scaling theory of turbulence and a comparison between the statistical approach to 
turbulence and field theory, we discuss a simple model of turbulent advection which 
is exactly soluble but whose exact solution is still difficult to analyze. The model 
exhibits a restricted universality. Its correlation functions contain terms with uni- 
versal but anomalous scaling but with non-universal amplitudes typically diverging 
with the growing size of the system. Strict universality applies only after such terms 
have been removed leaving renormalized correlators with normal scaling. We expect 
that the necessity of such an infrared renormalization is a characteristic feature of 
universality in turbulence. 



1 Introduction 



One of the basic and most successful ideas in theoretical physics has been that of univer- 
sality. It states that many systems with large or infinite number of degrees of freedom 
in certain asymptotic regimes display similar behaviors falling into general types. Such a 
situation has been encountered 

in statistical mechanics of 2 nd order phase transitions where the universality applies 
to the long distance asymptotics of correlation functions characterized e.g. by 
scaling exponents, independent of numerous microscopic details of the systems ]I| , 

in field theory where the universality implies the cutoff independence of the effective 
low energy description [[!]], 

1 Lectures read by the first author at the 34 th School of Theoretical Physics in Schladming, Austria, 
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in the theory of 1 st order phase transitions exhibiting universal shapes of critical 
droplets 0, 

in many-body condensed matter where universality limits possible behaviors of 
superconductors || or Hall fluids pj, 

in dynamical systems where it describes types of behavior of families of maps under 
iteration [^], 

in non-linear PDE's where it determines space and/or time asymptotics of solutions 
!• 

The basic idea behind universality of certain behaviors of systems with many degrees of 
freedom is that such behaviors are governed by few relevant degrees of freedom with simple 
dynamics. Those may be often exhibited using renormalization group transformations 
wiping out the irrelevant details of the system @ 1| . 

The question addressed in the present lectures is to what extent such a picture ap- 
plies to the fully developed hydrodynamic turbulence. It should be mentioned at the 
very beginning that the problem is still wide open although the first theory arguing for 
universality in turbulence, due to Kolmogorov ||, dates from 1941. In the 1 st lecture, 
we shall sketch how the main claims of the Kolmogorov theory relate to the cascade pic- 
ture of energy transfer in a turbulent flow. In the 2 nd lecture, we shall briefly present 
the functional-integral approach to the statistical theory of turbulence governed by the 
Navier-Stokes equation trying to stress the similarities and the differences with field the- 
ory. In the 3 rd and the 4 th lecture, we shall discuss a simple model of turbulent advection 
of a scalar quantity, known under the name of passive scalar (PS). The model is a good 
playground for testing the idea of universality in turbulence since one may obtain closed 
expressions for the stochastic initial data evolution and for the steady state correlation 
functions of the scalar. In the absence of external sources, the PS with deterministic ini- 
tial data undergoes a non-universal diffusion which averaged over initial data distributed 
in a homogeneous way becomes a universal super-diffusion. In the presence of random 
sources, the steady state 2-point correlator exhibits an energy cascade of the type expected 
in the Navier-Stokes turbulence. The exact solution for the higher point correlators are 
more difficult to analyze. We examine the origin of possible violations of strict univer- 
sality of these correlators. What emerges is a picture of restricted universality applying 
to infrared renormalized correlators from which few universal terms multiplied by non- 
universal coefficients where removed. In Conclusions, we summarize our discussion and 
point to possible directions of the future research. 



The PS model, although introduced long time ago jLO), has become recently a subject 
of intensive study, see |l I j |12| and also |13| (the latter, for a mathematical work on a 
simplified version of the model). In the final stages of the work on these notes, we have 



received preprints [[T4]][|T5| and ||16|| . The analysis of the first two has large overlaps with 



ours although we disagree with some conclusions of [14| and differ in some interpretations 



from 15 
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2 Energy cascade and Kolmogorov scaling 



The evolution of the local velocity field v(t,x) of the incompressible fluid at points x of 
the 3-dimensional space is described by the Navier-Stokes (NS) equation 

<%v + (v- V)v-z/Av = -Vp + f (2.1) 

supplemented with the incompressibility condition 

V- v = . (2.2) 

The presence of the dissipative term multiplied by the viscosity v of the fluid (with the 
dimension le ™ gth ) distinguishes the NS equation from the Euler equation. f(i, x) denotes 
the external force. p(t, x) is the pressure (divided by the constant density) and may be 
eliminated from both equations: 

Pit, x) = JL J v^H^l ^ . (2 . 3) 

If the flow takes place in a finite volume, boundary conditions should be specified. The 
size of the volume introduces an "integral scale" L into the problem and we shall assume 
that the external force f acts only at length scales comparable with L. Alternatively, we 
may consider flows in infinite space excited by external force with Fourier components 
concentrated in wavenumbers k with k = |k| < L^ 1 . Under the rescaling 

v(t, x) = erv(r£, sx) , 

f(t,x) = arf(rt,sx) , (2.4) 

p(£,x) = (JTS p(Tt, SX.) 

with as/r = 1, the viscosity v is replaced in Eq. (|2.1|) by v = ts~ 2 v . It is then 
convenient to introduce a dimensionless quantity, the Reynolds number 

R=¥, (2.5) 



where V stands for the typical value of the velocity differences over the integral scale L. 
R is invariant under the rescalings (|2.4j ). 

The basic phenomenological observation is that the hydrodynamic flows have very 
different behavior for small values of R (of order 1, say) where flows are regular (laminar) 
and large values of R (of order 10 3 and more) where the flows are very chaotic (turbulent), 
with complicated phenomena occurring at intermediate values of R. Here, we shall be 
interested in the regime of fully developed turbulence with very large Reynolds numbers, 
ideally in the limit R — > oo . It is sensible to use a statistical description of the flows 
in this regime. The stochasticity may be introduced into the description by considering 
random external force f or studying the evolution of random initial data or both. In a 
statistical description, the basic objects to look at will be the velocity correlation functions 
given by the expectation values of products of components of velocities. We shall mostly 
look at equal time correlators 



< 
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n=l 



"(t,X n ) > 



(2.6) 



which, in the stationary state of turbulence sustained by steady external forces, should 
be time independent. Besides, we may expect that, far from the fluid boundaries, the 
correlators (|2.6|) are translationally and rotationally invariant, the property that expresses 
the homogeneity and isotropy of the fully developed turbulence. The most commonly 
studied velocity correlators are the so called structure functions 

S N (x) =<(v(x)-v(0)) Af > . (2.7) 



One of the simplest consequences of the NS equation (|2.1| ) may be obtained by taking 
its scalar product with v and integrating the result over the space. These results in the 
relation 

|l/v 2 = ~/(Vv) 2 + /vf (2.8) 



which expresses the energy balance: the time derivative of energy on the right hand side 
is equal to the difference of the injection rate ff ■ v and the dissipation rate f /(Vv) 2 
of energy. Taking averages in the stationary state, we obtain 

<vf > = <^(Vv) 2 > = e, (2.9) 



i. e. the equality of the (intensive) mean injection and the mean dissipation rates of energy. 
How is this energetic balance realized over different modes of the fluid motion? The energy 
injection takes place at the distances of order of the integral scale by induction of big scale 
L eddies. According to the picture of the turbulent flow proposed in 1922 by Richardson 
17 j, the big eddies induce smaller eddies which, in turn, induce still smaller eddies and 



so on transferring energy from large to small distance scales. This process should not 
lead to a loss of energy until sufficiently small distance scales, say, smaller than rj, are 
reached. On scales smaller than 7] , the dissipative term v Av of the NS equation becomes 
important. One may formulate a more quantitative version of this picture. We shall follow 
here the discussion of the treatise fl8| which is an excellent reference to the problems of 



well developed turbulence. The idea is to localize the velocity field in the Fourier space. 
For v(k) = / e ~ ik ' x v(x) d 3 x define 



v <K (x)= / e tkx v(k)-^. (2.10) 

\k\<K K ' 



A more detailed version of the energy balance equation ( |2.8|) reads 

TtU< = -^/(Vv,J 2 + /v^ • f - JU K , (2.11) 



where U K is the density of energy flux out of the wave numbers with |k| < K . An explicit 
expression for 11^ may be obtained from the NS equation (2A) and may be argued to 
involve mainly v(k) with |k| = K , i.e. to be aproximately local in the wavenumber 
space. Taking the averages in the stationary state, we obtain the relation 



<v <K f > = <^/(w < j 2 > + <n A .> 



(2.12) 



which states that the mean injection rate of energy into the modes with |k| < K is equal 
to the mean rate of dissipation of energy e <K = < | /(Vv <Jf ) 2 > in those modes plus 
the mean flux of energy out of modes with |k| < K . According to the cascade picture 
of the energy transfer, the dissipation rate e K is essentially zero for K <C r]~ l . Since 
the injection takes place only around the integral scale, the change of the mean injection 



rate <p <K = < v <K ■ f > should be negligeable for K ^> L 1 so that <p <K should be 
equal to its K = oo value e for such K . It follows from Eq. ( |2.12| ) that in the so called 
inertial range of scales L^ 1 <C K <C r/ -1 , the mean energy flux in the wavenumber 
space ii k = < II K > is constant. This is an important quantitative consequence of the 
cascade picture of developed turbulence which, however, is not easily verifiable directly. 

Kolmogorov went further, postulating that in the inertial range ?) < x C I, the 
velocity structure functions may be expressed solely as universal functions of the mean 
dissipation rate e and the distance x , due to the locality and self-similarity of the cascade. 
The shape of the structure functions is then dictated by the dimensional considerations. 
One obtains 

S N (x) = C N e N/3 x N/3 (2.13) 

with universal coefficients CV (the right hand side is the only function of e and x with 
dimension ( ^" gt ) Ar ). Indeed, this behavior is consistent with the cascade picture. With 
the typical velocity of size / eddies vi ~ Sn(IY^ n ~ e 1 / 3 / 1 / 3 , their typical energy density 
jvf ~ e 2 / 3 / 2 / 3 and their typical turnover time t\ ~ l/vi ~ e -1 / 3 / 2 / 3 , one may estimate the 
energy flux from scales longer than I to those shorter than I as ~ vf/ti ~ e with no scale 
dependence. 



The scaling ( |2.13 ) of the 2 nd structure function translates under the Fourier transform 
to the relation 

E(k) =< v(k) • v(-k) > = / e~ ik ' x S 2 (x) rf 3 x ~ e 2/3 k~ u/3 (2.14) 



4irk 2 



so that the Kolmogorov theory prediction for the energy spectrum E(k) is 

E{k) ~ e 2 ' 3 k-^ 3 . (2.15) 

The theory allows also to locate the scale rj at which the dissipative effects become 
important, called usually the Kolmogorov scale. Estimating the dissipation rate on scale 

I as ~ v{— ) 2 , 77 should correspond to the scale at which v(— ) 2 = e. Using the relations 
I V 

v v ~ e^r/ 1 / 3 , v L ~ t^L 1 ' 3 , (2.16) 

we infer that 

V - = R'^L . (2.17) 



As should be expected, the Kolmogorov scale decreases with growing Reynolds numbers, 
i.e. with decreasing viscosity. In the limit R — > oo , the inertial range should invade 
all short scales. It does not mean, however, that the behavior in that range may be 
described by the Euler equation: the presence of even very small viscosity is essential 
for the existence of the stationary state of the turbulence by providing a mechanism for 
removal from the system of the energy injected by external forces. 

The verifications of the Kolmogorov theory is made difficult by the fact that both 
in experiments in the atmosphere, in wind tunnels or in liquid flows and in numerical 
simulations it is hard to obtain inertial ranges extending over many orders of magnitude 
necessary to extract the characteristic exponents of the structure functions. Neverthe- 
less, it seems that the behavior of low structure functions well agrees with Kolmogorov's 



predictions whereas for higher structure functions with iV ~ 10 one observes values of 
exponents lower than predicted (by ~ 10% or more). A multitude of more or less ad hoc 
mechanisms has been proposed to explain the possible departures from the Kolmogorov 
scaling for the higher structure functions. The underlying idea is that of intermittency 
limiting the turbulent activity to a subset of temporal and spatial modes. We shall not 
discuss the intermittent models here referring the interested reader to 



3 Functional methods in turbulence 

The Kolmogorov theory has a phenomenological character. It is strongly related to the 
cascade picture of turbulence which was postulated essentially independently of the NS 
equations. It may be compared to the mean-field approach to critical phenomena which 
leads to correct results only in limited situations and which, in general, has to be replaced 
by a more sophisticated theory taking into account the role of statistical fluctuations. It 
is then essential to try to build a theory of fully developed turbulence starting from the 
basic equations. 

For any functional evolution equation of the type 

d t <$> = + F , (3.1) 

where F is a functional Gaussian process, a Martin-Siggia-Rose (MSR) formalism 
permits to obtain formal functional integral expressions for the expectation values of 
functionals A(<&) in the stationary state. This is done as follows. We rewrite 

<A($) >= jf JA($) 5(d t $ + - F) det(d t + S -^§±) e-^^^MDF, (3.2) 

where C is the covariance operator of the Gaussian process F and M is the normalizing 
factor given by a similar integral with A — 1 . The role of the $ integral is to express 
$ as a solution of the evolution equation ( |3 . 1| ) (with a fixed initial condition). If the 
functional T(<&) is local in time then, formally, 

det(flfc + ^) = det(ft) det(l + df 1 ^) = det(ft) a**"™"* (3.3) 

since <9 t _1 (;£,£') = 6{t — if) and 5jF($)/5$ is proportional to 5(t — t') . Alternatively, 
one may use the (Berezin) functional integral over anticommuting fields to exponentiate 
the determinant. Rewriting also the delta function <5(<9 t $ + ^ r ( ( &) — F) as an oscillatory 
integral, we obtain 

<A($)> = -L j A ^) e i{A,d t ^n^-F) + ^ r 5H^-W^) DAD^DF . (3.4) 
Finally, the Gaussian integral over F leads to the relation 

<A($) >= ± J A{$) e i^,a t * + H*)) + htr&H*)/^-\(A,cA) DAD $ (3 _ 5) 

The MSR formalism applies to the Langevin equation describing the approach to 
equilibrium in statistical mechanics and Euclidean field theory. In this case, the functional 
JF($) appearing in the evolution equation ( |3.1| ) is of the gradient type: 



(3.6) 



and the Gaussian process F is the white noise, i.e. its covariance C is the identity 
operator. For example, for the scalar 4 theory, 



S{$) = J (f V$) 2 + \m 2 & + A$ 4 ) dt rf d x . (3.7) 



Perturbative treatment of the functional integral ( |3.5p based on expanding the exponential 
of the terms of order higher then 2 in $ and A into a power series gives rise to the 
perturbative expansion for the correlation functions. The latter reduces for the equal- 
time correlators to the perturbative expansion for the equilibrium expectations in the 
Gibbs measure e s ^ Dcf) with <\> as $ but without the time dependence. 



One may rewrite the NS equations fl2.1|) and fl2.2|) as a single evolution equation 

d t w = -P(v ■ V)v + z/Av + Pf (3.8) 

in the space of divergence-free vector fields, where P stands for the orthogonal projection 
on such fields. Taking the external force f random Gaussian, we end up in the setup 
formally resembling the Langevin dynamics. The important differences between the latter 
and Eq. (|3.8| ) are, however, that, first, the term P(v- V)v is not of the gradient type and, 
second, that we want to look at the noise f with the Fourier components concentrated in 
low wavenumbers, i.e. with the kernel of the covariance in the position space close to a 
constant as opposed to the delta function kernel for the Langevin equation. As a result, 
the role of the noise in the NS equation is to inject at small wavenumbers the energy which 
is then transferred to larger wavenumbers by an essentially deterministic process leading 
to the non-zero energy flux in the stationary state whereas in the Langevin dynamics, the 
noise plays an essential role on all scales resulting in the thermal state with no non-zero 
energy fluxes between scales. 



It is sometimes pointed |20[ that one encounters non-vanishing momentum-space fluxes 
also in the "equilibrium" field theory in the presence of anomalies. For the sake of illustra- 
tion, let us consider the chiral anomaly in two space-time dimensions. The chiral charge 
density of the Dirac field is 

J 50 (t,x) = ^xhWM) • (3.9) 

The anomaly calculation shows that in the external abelian gauge field in two space-time 
dimensions 

| <JJ 50 (t,x)dx>< K = ±fe^F„ u (t,x)dx, (3.10) 

and the right hand side is K independent. Here the subscript < K denotes a gauge- 
invariant (e.g. Pauli-Villars) subtraction at momentum- scale K . In other words, the flux 
of the chiral charge out of the modes with momenta < K is constant in K . It is not 
clear what lessons may be drawn from this analogy for understanding the inertial range 
of the fully developed turbulence. 



Applying the MSR formalism to Eq. ( |3.8j ), we obtain for the velocity correlation func- 
tions the expression 



N „ N 



< f[v^(t n ,X n ) >=j?J n^,X n ) eiMv+i^-.A^^Ca) jDaDV) (3n) 
n=l n=l 

where the a- and v-integrals are over divergence-free vector fields ( tr 5 drops out as 
it vanishes because of symmetry reasons). Separating the expression under the exponential 



into the quadratic part i(a, dtv — v Av) — |(a, Ca) and the cubic one i(a, (v • V)v) 
and expanding the exponential of the latter into the power series, we may generate the 
perturbation expansion for the velocity correlators, just like for the correlation functions 
of the Langevin evolution. The terms of the expansion may be represented by Feynman 
diagrams with the propagators read off the quadratic part under the exponential and the 
vertices off the cubic part. A closer look into the perturbation expansion generated this 
way shows that it is plagued by infrared and ultraviolet divergences (the latter in the limit 
of vanishing v ) . Various resummation techniques, mostly inspired by approaches used in 
field theory, were tried with the aim to improve the convergence, without complete success. 
Let us mention here another approach, originally labeled quasi-Lagrangian, which seems 
more promising. This approach has been developed over years by the Novosibirsk and 
more recently the Weizmann Institute schools f2"lf ||22|| [P3l ■ 



The idea consists of describing the turbulent flow in the frame moving with one of its 
points. Suppose that, for given velocity field v(t, x) , x(t) is the solution of the ODE 

f = v(f,x) (3.12) 

with the initial condition x(io) = xo . x(i) is the trajectory of the material point of the 
fluid located at x at time t . The velocity field, force and pressure in the frame moving 
with the fluid point are 

v'(t,x) = v(t,x + x(i)) , (3.13) 
f(t,x) = f(f,x + x(t)) , (3.14) 
p'(t,x) = p(t,x + x(t)) , (3.15) 

and they satisfy the equations following from the NS equations for v : 

d t V + ((v' - v' ) ■ V)v' - uAV = -Vp' + f , (3.16) 
V-v' = 0, (3.17) 



where v' (t) = Vq(£, 0) . Eqs. (|3.16|) and (|3.17Q may be again written as a single evolution 



equation in the space of divergence-free vector fields: 

d t V = -P((v'-v' ) ■ V)v' + uAV + Pf (3.18) 

to which the MSR formalism may be applied if we assume that f is a Gaussian random 
field. This is not the same as a similar assumption about f but as long as the Fourier 
components of the force are concentrated in small wavenumbers, the details of the force 
distribution should not matter in the inertial range. In such a situation, one may expect 
that the equal-time velocity correlators in the stationary state of the evolution equation 
( |3.18| ) coincide in the inertial range with those corresponding to the equation fl3.8| ). The 
perturbation expansion generated by the MSR technique for the v'-correlators has the 
same propagators but modified vertices which break the translation invariance and make 
this way the quasi-Lagrangian expansion more difficult to analyze. Nevertheless it seems 
that the quasi-Lagrangian approach leads to drastic improvement in convergence of the 



(Schwinger-Dyson resummation) of the perturbative expansion [21] [23] and it is possible 
that it provides an important clue to the understanding of the statistical properties of the 
turbulent flow. 



4 Passive scalar 



4.1 Definition of the model 



In view of the difficulties encountered by the theory of the Navier-Stokes turbulence, it 
may be useful to search for simpler systems exhibiting some of the aspects of full-fledged 
turbulence but easier to control. Such a system seems to be provided by a model of 
passive advection in a random velocity field v(t, x) of a scalar quantity T whose density 
T(t, x) satisfies the equation 

d t T + (v-V)T - uAT = f , (4.1) 

where now v denotes the molecular diffusivity of the scalar T and f(t, x) describes the 
external sources. In the ideal situation, we should consider v as the velocity of the NS 
turbulent flow, but we shall, instead, assume that v(t, x) is a centered Gaussian field 
with the covariance 

< v* (t, x) v 3 (f, x') > = 6{t - D ij (x - x') , (4.2) 

i. e. white noise in time. The spatial part of the covariance will be taken, for concreteness, 
as 

D y (x) = D J e ik x (k 2 + m 2 )-( 3+K )/ 2 (5 ij - k'k j /k 2 ) ^ , (4.3) 

where the transverse projector in the Fourier space ensures the incompressibility of v. 
Small m 2 should be viewed as an infrared cutoff making the integral convergent for < 
k < 2 . Note that the dimension of Dn is len ? th — . For the 2 nd velocity structure function 

w time J 

we obtain 

< (v(t, x) - v(t, 0)) 2 > ~ D u (0) - D*(x) = D 4i (x) (4.4) 

which, for < k < 2, has the m — > limit proportional to |x| K . i.e. growing with the 
distance. More exactly, 

Dij w = iS^k 611 D ° m ~ K » ( 4 - 5 ) 

i. e. it diverges with m — > , but 

lim Zfr'(x) = Di f(2 + /t)<5 ii |x| K - /txV|x| K - 2 ) (4.6) 

and it is a homogeneous function of x . D\ = 22 +i^ 2 372^[ 3 2 | K ) D . Note the singularities at 
and 2 in the K-dependence. The Kolmogorov scaling ( 2.13|) corresponds to k = | which 



lies in the interval < k < 2 that we shall consider. It should be mentioned, however, 
that the time-decorrelation, evenness and Gaussian character of the velocity distribution 
is a rather drastic departure from the statistical properties of the turbulent NS velocities. 
Nevertheless, the model, dating back to 1949 pTOfl , possesses an interest of its own and 



has been used to describe a variety of phenomena like poluant or tracer transport and is 
closely related to models describing forced diffusion through porous media ||24|| . 



We shall take the field /(t, x) describing the external sources also random Gaussian, 
independent of v, with mean zero and covariance 

< /(*, x) f(t', x') >=S(t- C(^) = S(t - C L (x - x') (4.7) 



where C is a real positive-definite function from the Schwartz space <S(R 3 ) and L is the 
integral scale. Sometimes, it will be useful to consider also the advection of a complex 
scalar in which case C may be a complex positive-definite function. We would like to 
study the statistical properties of the solutions of Eq. Q4.1|) in the regime of small u, 
small m (which may viewed as the inverse of another integral scale) and large L (in the 
limit L —>■ oo, the field / has only the zero wavenumber component). In particular, the 
universality question for the PS may be formulated as inquiring about the existence of 
the limit of the correlation functions 

N 

< n^n,x n ) > 

n=l 

in a stationary state of the system when u, m, L" 1 — > and about the independence of 
such a limit of the shape of the source covariance C . 



(4.8) 



4.2 Exact solution 

The big advantage of the PS model over the NS case is that one may obtain closed 
expressions for the correlation functions of the scalar (the assumed time-decorrelation of 
the velocity field plays here the crucial role). For 2-point function, this was first done in 
[p5fl 3 . In this sense the model is exactly soluble. In order to work out its exact solution, 
let us rewrite the basic evolution equation ( [4.1| ) of the scalar in a general form 

d t T + (3T + vaT = f , (4.9) 

where (3 is a skew-symmetric random operator ((3 = v • V in our case) and a is a 
positive operator (a = — A). The (3 term in the equation is conservative whereas the a 
one is dissipative (for the energy \ ||T|| 2 ) . For the sake of simplicity, we shall argue in the 
finite-dimensional case, i.e. when T(t) and f(t) take values in R d and f3(t) and a are, 
respectively, skew-symmetric and positive symmetric d x d matrices. For the covariances 
of independent centered Gaussian processes (3(t) and f(t) , we shall write 

< (3ij(t) Ph(s) > = 5(t-s) B ijM , 

<fi(t)f j (s)> =6(t-s)C tJ . (4.10) 

In finite dimensions, the original stochastic PDE becomes a stochastic ODE which is 
easier to handle. Let us simplify the matter further by assuming, for a moment, that 
(3(t) and f(t) are continuous functions. Then the solution of Eq. ( |4.9|) with the initial 
condition To at t = to takes the well known form 

T{t) = R(t,t )T + [ t R(t,s)f(s)ds , 

where R(t,to) is given by the ordered exponential (t >t ) 

R(t,to) = Te-^ [va+m)dT (4.12) 

(T stands for the time ordering). R(t, to) may be also expressed by the following limiting 
formula : 

R(t,t ) = lim exp[— / [ya + (3{r)) dr\ exp[— / [ya + (3{t)) dr] ■■■ 



3 We thank U. Frisch for explaining to us the 4-point function case 



(4.11) 



••• exp[- f(ua + /3{T))dr] = lim R Q (t,t Q ) , (4.13) 

where (I q ) is a time ordered partition of the interval [to,t] into subintervals of length 
(t-t )/Q. 

In the random case, the white noise stochastic processes f3(t) and f(t) do not have 
continuous realizations. They need smearing to give genuine random variables. It is easy 
to see that / J( 5(r)(ir = (3 1 where / are intervals of time, are already genuine Gaus- 
sian random variables which are mutually independent for non-intersecting or adjacent 
intervals: 

<(3 I l] (3 J kl >= \InJ\B ij>kl . (4.14) 
Let us denote by B the covariance operator acting on skew-symmetric matrices by 

(B(3)ij = T,B ijM Pki ■ (4.15) 

ijkl 

Necessarily, B > 0. We shall also need below a contracted version B of B with the 
matrix elements 

B *i = 4E^W • ( 4 - 16 ) 

k 

B is a symmetric operator acting in R d . Note that 

< ((P 1 ) 1 (5% > = -< (/3 7 /3% > = 2|/|B y (4.17) 

from which it follows that B > 0. In the random case with white noise /?(£), the 
approximate evolution operators RQ(t,t Q ) under the limit on the right hand side of Eq. 



( [4.13| ) still make sense for each Q as random (matrix valued) variables and one can show 



that 
Proposition 1. 



lim < R Q (t,t ) > = e -(*-*o)(^+B) = e -(*-*o)Mi . (4.i 8 ) 

Q— >oo 



In order to indicate how the above result arises, consider the perturbation expansion for 
R(t, t ) in powers of (3 



R{t,t ) = El- 1 )™/ dr m j dr ro _!--. j dn ^- Tm)va I3{r m ) e~^~^> a /3(r m ^ 

■■■ P(n) e - (Tl - i(,) ^ . (4.19) 



m=0 to t t 



In the white noise expectation value of the right hand side computed with the use of the 
Wick Theorem, only neighboring (3(r m ) may be paired forcing the respective times r m 
to be equal. One obtains this way the perturbation expansion of e _ (*~*°)^ Q!+ ' B ) in powers 
of B . Proposition 1 is a rigorous version of this result, not very difficult to prove p6| . 



Similarly, one may calculate expectations of products of matrix elements of Rq(t, t ) . 
We shall use the tensor product notation Rq(t, t )® N as a bookkeeping device for all such 



products. Let us denote by /C the symmetric operator on R^C^R^ — Matdxd built of the 
coefficients so that 

{Kg)ik = ^Bij,u9ji ■ ( 4 -20) 
H 

Note the difference between JC and the covariance operator B of the process /3(t) defined 
in Eq. ( 4.15| ). In what follows an important role will be played by the iV-body symmetric 
operators M. N acting in (R d )® N = K Nd defined by 

N 

M N = (M x )n ~ E (^W , (4.21) 

n=l l<n<n'<N 

where M.\ = va + B , {M.\) n denotes M.\ acting in the n th factor of the tensor product 
(R^)®^ and, similarly, {K) n , n t stands for K, acting in the n th and n' th factors of (R^)®^ 



Proposition 2. For N = 2, 3, . . . , 

lim < R Q (t,t )® N > = e - (t - to)A ^ . (4.22) 

Q^oo 

Again, this result is not surprising if we use the perturbation expansion (|4.19|) for each 
R(t, to) and note that the white noise expectation may pair now, besides the nearest 
neighbor f3{j m ) in the perturbative terms of the same R(t,t ) also /?(r m ) corresponding 
to different R(t,to), resulting in the perturbative expansion of e~^ t ~ to ' MN in powers of 
operators B and /C. In fact, one may show [^6| that -Rq(t,to) converge to a limit in 



all L p spaces, p < oo , of the white noise (5 . Naturally, we shall denote the limit by 
R(t, t ) . R(t, to) are contracting operators, they are continuous in t, to in the L p -norms 
and satisfy the composition law. We shall take then T(t) given by Eq. ( 4.1 1| ) with the 



limiting R(t, to) as the solution of the stochastic equation (|4.9|) with initial value TJ 



o ■ 



In the absence of the external sources the time evolution of the equal-time correlators 
of the scalar T simply reduces to the relations 

< T(t)® N > = e-<*- t °) M »T(t ) 9N , (4.23) 

hence it is given by the semigroups (e~ tMN )t>o which may be studied by looking at the 
Fourier transforms of the semigroups i.e. at the resolvents (iu + M.^)' 1 ■ In any case, 
the properties of the time evolution of the statistics of initial data are determined by the 
spectral properties of the symmetric operators Mn ■ It is easy to see that the operators 
M.N satisfy the inequality 

N 

M N > vJ2( a )n ■ (4-24) 

n=l 

An important consequence of this inequality is the exponential time decay of the corre- 
lation functions < T(t)® N > as long as the dissipation is present in the equation ( |4.9| ) 
(recall that we are looking at the finite-dimensional case). The decay rate is given by the 
lowest eigenvalues Ajv,o of JAn and is bounded below by iVz//i where //o > is the 
lowest eigenvalue of a. If we remove the dissipative term from the equation by setting 
v = then the long time behavior of < T(t)® > will have a persistent tail since the 
operators M. # develop zero modes when v is taken to zero. Suppose now that the initial 



data T(to) are also random and independent of (3{t) for t > to- Averaging Eq. ( |4.23j ) 
over T , we obtain 

< T(t)® N > = e - (t -* o)A1 ^ < T(t )® N > . (4.25) 

A natural problem is the existence of an invariant measure on the space of initial data. It 
follows from relations ( |4.24| ) and ( |4.25| ) that there are no invariant measures with finite 



moments on the space of initial data as long as v > . On the other hand, when v = , i. e. 
in the absence of dissipation, one may show that any measure preserved by a. a. rotations 
e 13 is invariant and that this property characterizes the invariant measures with finite 
characteristic functions |26 . In particular, the Gibbs measures ~ e -const.||T|| ^dji are 



invariant under the flow (L9) when v — 0. As we shall see, the situation is quite different 



for infinite number of degrees of freedom. 

In the presence of white-noise sources and of dissipation, the 2-point equal time cor- 
relation function evolves according to the equation 

< T(t)® 2 > = e - {t - to)M2 T{t )® 2 + jds e - {t ~ s)M2 C 

to 

= e- {t - to)M2 T(t f 2 + (l-e'^-^^M^C (4.26) 

which is a solution of the differential equation 

dt < T(t)® 2 > = -M 2 < T(t)® 2 > + C . (4.27) 

When t — > — oo, the term with T(t ) disappears due to the positivity of M.2 (for v > 0) 
and we obtain 

< T m > = M^C . (4.28) 

This is the time independent 2-point function in a steady state in which the quantities 
of the energy injected by the force and dissipated balance each other. Similarly, for the 
general equal time correlators in the presence of the external sources and of dissipation, 
the dependence on the initial conditions is wiped out when to — ► — oo due to the positivity 
of all M. n ■ The higher stationary state equal-time correlation functions take the form: 

< j-,®2M+l > _ q 

<T® 2M > = aF ™, (4.29) 

where F^m £ (j^)<x> 2M j s a single 2-particle channel contribution. Permutations o in- 
terchange the factors in (R d )® 2M , Sm is viewed as a set of permutations of the pairs 
(2m — 1, 2m) , m = 1, . . . , M , and S^ 1 as that of permutations acting within the pairs 
so that the sum is over the pairings of {1, . . . , 2M} or 2-particle channels. 

F 2M = Mt£t E ° ((>W-2 ® 1 2 ) • • ■ (M, 1 ® Um-2)) C 0M . (4.30) 



a-es 



M 



Exact expressions for non-equal times A^-point correlators are also easy to obtain. They 
involve also the heat operators e~ tMN ' with N' < N. 



Let us discuss briefly some multibody combinatorics behind the solution ( [4.30|) . Let 



us denote by Vim the set of pairs (2m — 1, 2m) , m = 1, . . . , M . We shall define the 
connected parts F 2 C M of the single channel correlator by the inductive formula 

F2M = E 0^ ( 4 - 31 ) 



partitions 7r£ll 

n of Tom 



in the, hopefully, self-explanatory notation. The connected 2M-correlator < T® 2M > c is 
given by Eq. (|4.29|) wi th F 2 C A/ replacing F 2M on the right hand side. Note that, by virtue 
of the definition (|4.21| ) of M.n , for each partition II of V 2 m , we may write 

M 2M = E ( 4 - 32 ) 

?ren {7r,7r'}cn 

where A4 n is A^2|tt| acting on the tensor factors of (R, d )® 2M labeled by the pairs of tt 
and 

Cn,*' = E (fyn,n' ■ (4.33) 



The basic equation for the connected single 2-particle channel correlation functions K 
is given by 



c 

2M 



Proposition 3. For M > 1 , 

M 2M F 2 C M = E (4.34) 

n={7r,7r'} 



Note the connected character of the right hand side. Let us prove this relation. From the 
definition ( [4.30| ) of the single channel function, 



M 



2M r 2M 



E c -® F -' ■ 



(4.35) 



n={7r,7r'} 
|ir| = l 



■M-2M F 2M 



On the other hand, using Eqs. (|4.31| ) and ( [4.32| ), we obtain 

= M 2M F 2M - J2 K 

n TrGn 

|n|>i 

= M 2M F 2M - E 52( M «K)( 
n Tien 

n|>i 

+ E E (^Af;®f;,))§. 

n {7T,7r'}cn 
|n|>i L ' 



ir'en 



(4.36) 



7r"^7r,7r' 



Notice that on the right hand side, in virtue of Eqs. (|4.35|) and ( |4.28|) , the 1 st term cancels 
with the contribution to the 2 nd one from 7r consisting of a single pair. Rewriting M. n F^ 
for the [-7T | > 1 contributions to the 2 nd term using Eq. ( 4.34 ) as the inductive hypothesis, 
we observe that all these contribution are cancelled by the 3 rd term. What is left from 
the 3 rd term after the cancellations gives exactly the right hand side of Eq. ( P~B~ 



Although we worked directly with the evolution equation ( f4.9[ ) in order to obtain 
expressions (|4 . 2 8|) and (|4.30| ), one may obtain them as well using the MSR formalism 
which gives 



< rp®2M > 



1_ 

77 



T 



)2M e i(A,(dt+va+f3)T)-^(A,CA)-i 



(P,^ 1 (3) 



DA DT DP . (4.37) 



Performing the Gaussian A- and T-integrals, we obtain 



<T 82M >=i ]T I l[(d t + va + (3)- 1 C(-d t + va- (3)- 1 e'W 3 ) Df3. (4.38) 

parings pairs 

of {1,...,2M} 

If we expand 

oo 

(dt + ua + P)- 1 = J^i-indt + va)- 1 (P{dt + va)- x ) m (4.39) 



m=0 



and use the fact that the kernel of (dt + va)' 1 is equal to 9(t — t') Q-'y t - t ') ua ) the /3-integral 
in ( [4.38D reduces to the previous calculation and it reproduces the expressions (|4.28j) for 
the 2-point function or Q4.2S ) and (|4.30| ) for the 2M-correlator. 



4.3 Eddy diffusion versus super-diffusion 

The exact solution for the PS correlation functions obtained above in the case of finitely 
many degrees of freedom still makes sense for the infinite-dimensional case of the evolution 
equation (|4.1| ). Instead of starting with this equation, and extending the previous analysis 
from the stochastic ODE case to the stochastic PDE one, which would require additional 
work, we shall directly study the infinite-dimensional version of the exact solution (|4.29| ) . 
( f4.30|) involving operators rather (generally unbounded) than finite matrices. The analysis 



of the expressions for the PS correlation functions obtained this way is more complicated 
than in finite dimensions, but also more interesting. 

The covariance B of the Gaussian process (5 = v • V and the operator B related to 
the expectation of /3 2 , see Eqs. fl4.16| ) and ( |4.17| ), are given by the kernels 



i3(x 1 ,yi;x 2 ,y 2 ) = V v (x x - x 2 ) ^(xi - yi) <9,-5(x 2 - y 2 ) , 

S(x,y) = ~\j £(x,z;z,y)d 3 z = -^'(0) c^(x - y) . (4.40) 

Since £>' J '(0) = ±8 ij V ll (0), B = -~V ll (0) A and, consequently, 

Mi = va + B = -{v + \V ll {$)) A = -u eS A . (4.41) 



It follows from Eq. ( |4.23j ) with N — 1 that, in the absence of the sources, the expectation 



value of the scalar diffuses 

< T(t,x) > = J e-( t - t °)- Ml (x,y)T(t ,y) d 3 y 



/e ( '-' o) ^ A (x,y)T(t ,y)rf 3 y (4.42) 



with the effective diffusion constant v c r composed of the molecular diffusivity v and the 
eddy diffusivity |X>"(0) . For small m , the eddy diffusivity dominates and the diffusion 
is driven by the large distance scales (recall from Eq. (|4.5| ) that V ll (0) = 0{Dqui^ k )) . 



Operator /C with kernel /C(xi, x 2 ; yi, y 2 ) = B(x.\, yi; x 2 , y 2 ) acts on functions g of 
six variables by 

(/C<7)(x 1)X2 ) = ^(x 1 -x 2 )9 1 . 1 ^^(x 1 ,x 2 ) . (4.43) 



so that for the iV-body operators M. n given by formula ( f4.21| ) we obtain 

N 

M N = -E(^ A x„ +I^'(0)^%) " E ^(Xn-V)^-^ 
n=l n<n' 
JV 

= -E ^ff A Xn - E ^'(xn-^jS^ . (4.44) 

n=l n<n' 

In particular, 

M 2 = _ z/(Al + a 2 ) - ip^OR,^ - i^(0)a 4 ^ - 2^( Xl - x 2 R^ . (4.45) 

Note that .M2 commutes with (three-dimensional) translations and in the action on 
translation-invariant functions of xi — x 2 = x reduces to 

M 2 = -2uA - V^i^didj . (4.46) 

Since X> y '(x) = P tJ '(0) — £> lJ (x) has an m -> limit, so does the operator .M2 in the 
action on translation-invariant functions and when v — > , it becomes a singular elliptic 
operator 



M S 2 C = -D 1 ((2 + K)8 ij \x\ K - Kx i x j \x\ K ~ 2 ) didj , (4.47) 
:((2-*)/2) 

2 2 + k 7t 3 / 2 k(3+k) 



where .Di = 22 +i| r 2 3 /2 K j^l ^ -Po ? see Eq. ( |4.6|) . The dimension of D\ , as that of D , is 



length 



equal to 

^ time 

Similarly for > 2 , in the action on the translationally invariant sector, 

TV 

m n = -v E A x„ + E &j (x n - x„o (4.48) 

n=l n<n' 

and has an m — > and v — > limit which is a singular elliptic operator . This 
has important consequences. Although the time evolution given by Eq. (Q5) of fixed 



initial conditions is dominated for small m by the eddy diffusivity diverging as m — > . 
the time evolution governed by Eq. (|4.25|) for random initial data with a translationally 
invariant distribution behaves well in the limit when m — > and v — > . In particular, it 
follows from the scaling properties of operators A4 S § that, if the translationally invariant 
initial moments < n.n^(^0;X n ) > are fast decaying in the difference variables (weakly 
correlated initial data) then, for m — , the rescaled correlations 

A2(ir7o < Yl T(\t, A~x n ) > (4.49) 

n=l 

have a non-trivial limit when A — > oo . In other words, the equal-time correlators decay 

__3N_ 

super-diffusively in the absence of sources, like t 2 ( 2 - k ) , i. e. faster than diffusively (like 

3JV _____ 

t~~ ). Similar behavior in a simplified model of passive advection was discussed in [fL3|]. 

Let us look closer at the operator .M 2 in the translationally invariant sector. Since 
the limit m — > presents no difficulty in that case, we may directly work in the m = 
case, i.e. with operator Ai 2 = —2uA + Aif , see Eq. ( |4.47| ). In the radial variables 

M 2 = -id r (ur 2 + D ir 2+K )d r + l ^{2v + D 1 (2 + k)^ 2 ) , (4.50) 



where I is the angular momentum {M.2 commutes with rotations). As for the resolvents 
{iuj + they may be explicitly computed in the translationally and rotationally 

invariant (I = 0) sector when v — > 0. In particular, one obtains for the kernel 



l + K „ . 1 , ,2 



(iw + A<f)- 1 (x,0) = uj— f{u*\x\ — ) (4.51) 

2(1 + k) 

with an exponentially decaying function /, /(y) ~ y 2 -" for small y. / may be 
expressed in terms of the Hankel functions. This result has an important implication. 
Recall that in finite dimensions the invariant measures for the stochastic evolution of the 
initial data were characterized by the symmetry generated by rotations e 13 . Suppose 
that, as for a finite number of degrees of freedom, the evolution ( |4.25| ) preserves in the 



limit v — > the white noise Gibbs measure ~ e const -ll r ll JJT with the 2-point function 

<T(x)T(y)>~<5(x-y). (4.52) 

This would imply that 

/ g(x) (iu + M?)-\x, 0) rf 3 x = - g(0) (4.53) 

for any test function g, which is in contradiction with Eq. ( |4.51| ). Taking Fourier trans- 
forms in space and time of the latter, one obtains 



- tkx 



e^^ c ( X;0 )d 3 x = p(^|k|) , (4.54) 



with p(0) = 1 and p decaying at infinity, instead of a constant as would be implied by 
Eq. (|4.53|) . One should compare the relation ( 4.54 ) to the relation 



-*k.x e -* A ( x , 0) d 3 x = e-* |k l 2 (4.55) 

describing the diffusive decay of the energy spectrum. For example, Eq. ( |4.54|) implies 
that for homogeneous, fast decaying initial data 2-point function, 

f <T(t,0)T(t,x) > x 2 d 3 x ~ (4.56) 

instead of being proportional to t : again a super-diffusive behavior. The reason of the 
failure of the Gibbs measure to be invariant under the evolution ( |4.25| ) for infinitely 
many degrees of freedom is the energy cascade towards degrees of freedom with smaller 
and smaller wavenumbers, see below. This failure shows that, for the infinite number of 
degrees of freedom, the symmetry group generated by the exponentials e 13 , which for 
(3 — v • V is the group of volume preserving diffeomorphisms of R 3 , is broken. It is not 
excluded, however, that this symmetry may still be useful in analysis of the model. 



4.4 Quasi-Lagrangian approach to passive scalar 

PS constitutes a perfect model to test the quasi-Lagrangian ideas mentioned at the end 
of Section 4. The scalar density T'(t,x) = T(t, x+ x(t)) viewed from a point moving in 
the fluid satisfies the equation 



d t T' + (v' - V )T' - 2/ AT' = /' 



(4.57) 



where v'(£, x) = v(£, x + x(i)) and /'(£, x) = f(t, x + x(i)) in the notation of Section 
4. It is again a stochastic equation of the type (|4.9| ) with j3 — (v' — Vq) • V. One 
may show that for Gaussian homogeneous v(t, x) and f(t,x), white noise in time, the 
distributions of the random fields v'(t, x) and /'(t, x) obtained by the quasi-Lagrangian 
transformation become the same as those of v(t, x) and /(£,x) asymptotically for large 
times. For v'(t,x) and f'(t,x) distributed as v(t,x) and /(i, x) before, we may produce 
the exact solution for the correlation functions of T" by replacing the operators M. ^ by 
new operators -M'^y obtained by changing the velocity correlation function (}4.2|) to 

< (t/'(t, x) - v'% 0)) (t/''(f , x') - v' j (t', 0)) > 

= 6{t - if) (l^'(x - x') - D ij '(x) - Z^'(x') + D ij '(0)) 

= - 5(t - /) (l) y '(x - x') - D ij (x) - D ij (*.')) • (4.58) 

Notice that only the subtracted functions with the regular m — > limit enter the 
correlator (|4.58|) . One obtains 

N N 
n=l n=l 

+ E P j (xn - x„0 - ^ W - ^ y '(x„0) ^ n 9^ . (4.59) 

n<n' 



The crucial observation is that in the action on translationally invariant functions oper- 
ators M.' N and M.n i the latter given by Eq. ( (4.48Q , coincide (the terms with D !3 (x„) 
cancel). Hence for the PS, the stationary equal time quasi-Lagrangian correlators will 
be equal to the original (Eulerian) ones for the homogeneous sources. This is not the 
case for the non-equal time correlators which involve the heat operators e~ tMN ( e~ tM N ) 
in the non-translationally invariant sector. The latter are governed by the eddy diffu- 
sivity, i.e. by the integral scale in the Eulerian case but show a universal behavior in 
the inertial range in the quasi-Lagrangian description. The quasi-Lagrangian picture is 
also better behaved from the point of view of the standard perturbative approach which 
consists of tre ating the 2-particle terms in M. n as a perturbation of the 1-particle ones, 
see Eq. ( 4.21 ). In the Eulerian approach, the unperturbed terms are dominated by the 
diverging eddy diffusivity whereas in the quasi-Lagrangian one the 1-particle terms have 
a regular limit when m — > 0, compare Eqs. ( |4.44|) and ( |4.59| ). Thus the quasi-Lagrangian 
perturbation expansion seems to be a better tool for the analysis of the behavior of the 
model. 



5 Steady state of the passive scalar 
5.1 Universality of the 2-point function 

Let us pass to the study of the stationary state 2-point function < T(xi) T(x2) > = 
F 2 (|xi— x 2 |) in the presence of the random sources with a rotationally invariant covariance 
Cl = C(- / L) , see Eq. (|4.7|) . F 2 satisfies the differential equation 



(M 2 F 2 )(r) = 



~d r (ur 2 + D ir 2+K )d r F 2 (r) = C L (r) 



(5.1) 



which is implied by Eq. ( |4.28| ). Integrating once, we obtain 



jC{p/L)p 2 dp 

= ~\ ° 2 I n 2+K > ( 5 - 2 ) 

where the constant of integration has been chosen so as to obtain d r F 2 vanishing at zero 
as required if F 2 is to be a function on R 3 regular at zero. One more integral gives 

oo/C(p7L)p ,2 c/p' 
J v p z + Dip z+K 

r 

where the choice of integration constant assures that F 2 decays at infinity. We want to 
study the behavior of F 2 (r) when the integral scale L — > oo. To this end, let us rewrite 
the p-integral from r to oo as a difference of the integral from to oo and from to r : 

y JC(p')p' 2 dp' r J C(p' / L) p' 2 dp' 

*W = \^ I + 5^ *> ~ I / VpTTD^ d " • (5 ' 4) 

r r r r 

The first of the integrals is a C-dependent constant diverging when L — > oo . On the 
other hand, the integral from to r has a limit when L — > oo equal to 

-»C(0)/ 2 f^, + , (5.5) 

which is universal (i.e. depends only on C(0) which, as we shall see in the next Section, 
equals twice the energy dissipation rate e of the scalar). Note that the non-universal 
term (a constant) is annihilated by M. 2 ■ This has to be so if Eq. ( |5.1D is to be satisfied. 
The right hand side becomes universal in the limit L — > oo so all non-universal terms 
in F 2 (r) surviving in this limit have to be annihilated by hA 2 . We shall see this general 
mechanism limiting possible non-universal terms also for the higher point functions. The 
constant term of F 2 is automatically subtracted in the 2 nd structure function: 



S 2 (r) = 2(F 2 (0) - F 2 (r)) = \e I , 7 (5-6) 



o 



p 3 dp 
up 2 + Dip 2 



where the last equality holds for infinite L. In the v — > limit, 

= r 2 - K (5.7) 

and the same universal result holds approximately in the whole inertial range ^<r<I 
where at the Kolmogorov scale r\ = 0((u / 'Di) l l K ) the term ur 2 becomes comparable to 
D\p 2+K . The 2 nd structure function exponent ( 2 = 2 — k. 

It is instructive to repeat the same analysis for anisotropic sources, i.e. without 
assuming that the correlation Cl is rotationally invariant. We shall have then to solve 
the equation 



M 2 F 2 = C L 



(5.8) 



in other angular momentum sectors. This is an easy exercise [2B] and one still obtains, 
for m and v equal to zero, the result 



- 7^ 2 - K ^ - i(^|x| 2 - K es F| c (x), (5.9) 

where 7 is a C-dependent constant. In the case of a complex PS, however, further 
subtractions are needed when < k < 1. For F 2 (x) = < T(0) T(x) >, one shows then 



2t that 



where real 7 and imaginary 7, are C-dependent constants. The 7j-terms are present 
only if C has angular momentum 1 components and this is impossible in the real case 
where only even I 's appear^. Again, the non-universal terms are homogeneous zero modes 
of A4 2 which annihilates constant and linear functions, in accordance with the general 
argument. Similarly as the constant, the linear terms appearing with the coefficient 
L 1 ^ K in the presence of anisotropic sources for the complex scalar drop out from the 2 nd 
structure function < |T(x) — T(0) | 2 >= 2F 2 (0) — F 2 (x) — F 2 (— x) and do not spoil the 
inertial range universality of the latter. 



5.2 Energy cascade 

Let us look at the energy balance of the scalar. The dynamical equation ( |4.27| ) takes the 
form 

d t F 2 (r) = -(M 2 F 2 )(r) + C L (r) . (5.11) 

and reduces in the stationary state to Eq. (|5.1| ). In particular, the mean energy density 
e = < \T 2 > satisfies 

d t e = ld t F 2 (0) = - 1 -(M 2 F 2 )(0) + ±C(0) 

= i/(AF 2 )(0) + iC(0) (5.12) 

(only the — 2z/A term of A4 2 contributes at r = 0). 

e = -u(AF 2 )(0) = -v-Ld r r 2 d r F 2 (r)\_ (5.13) 

is the mean dissipation rate and ^C(O) is the mean injection rate of energy (both inten- 
sive) . In the steady state the energy balance reduces to the equality of the rates : 

e = iC(0) . (5.14) 

It is, of course, satisfied by our solution for the steady state 2-point function. 

We would like to see how the energy is distributed among different wavenumbers 
following the schematic discussion of Section 1.2. For the mean energy density of modes 
with |k| < K , e<K = < %T<k >> we obtain 

e< K =\ I (/e-^F 2 (y)d 3 y) tSt = 1<k)dk (5.15) 
4 K.G. would like to thank R. Kraichnan for pointing this out to him 



with the energy spectrum 

e(Jfe) ~ k~ 3+K (5.16) 
in the steady state inertial range. The dynamical equation (|5.11|) implies that 

d t e< K = -e< K + (p<K - tt k , (5.17) 

where 

e<K = -v I (/e-^AF 2 (y)rf 3 y) tSt 

|k|<K V J ( n) 

= _2 u J sm(Kr)-Krcos(Kr) Q^^p^ dr (g^g) 
7T Q r 

is the mean energy dissipation rate in the modes with |k| < K and 

<p<K=\ I (f C L (y) d*y) = \ J zJSl C (r/L)dr (5.19) 

2 |k|<K V ' { > r 

is the mean energy injection rate into |k| < K . 

nR = J sHKr)-Krcos(Kr) drr 2 + . dr p^ r) dr (5 20) 



comes from the contributions of {M. S 2 C ^){ r ) to (|5.11|) . In view of the relation ( |5.17| ), 
it is natural to interpret ttk as the mean energy flux out of the modes with |k| < K , 
compare Eq. ( 2.1 1|) in Sect. 2. In particular, in the forced steady state, we obtain the 
energy balance equation 

<P<K = £<K + 7Tk , (5.21) 

compare Eq. ( [2.12|) . 

With the Fourier transform of C(- / L) concentrated around the zero wavenumbers, the 
energy injection is essentially limited to small wavenumbers and the mean injection rate 
ip<K is approximately constant for K 3> L^ 1 and equal to its K = oo value iC(0) . On 
the other hand, the mean dissipation rate €<k m &y be easily showed for k > 1 to satisfy 

e< K < 0(^-K K ) . (5.22) 
Indeed, it follows from Eq. ( |5.2| ) that \d r r 2 d r F2(r) \ < const. Df 1 r 2 ~ K . Estimating 

\sm(Kr)-Krcos(Kr)\ < h™st [Krf for r < K^ , ^ 

I const. Kr for r > K 



we obtain 



ve 1 / K oo \ 

e< K < const. — [K 3 J r 2 ~ K dr + K / r" K rfr (5.24) 

D\ \ l/K J 

which implies the desired bound ( |5.22|) for k > 1 . For k < 1 , a better estimate for 
\d r r 2 d r F 2 (r)\ for large r using the decay of C implies a somewhat worse bound 

6< x < O(g-K^) + O(^^K) . (5.25) 



As a consequence, for k > 1 , the mean dissipation rate e<K stays very small relative to 
e as long as K< r^ 1 where r\ = 0( (u/Di) 1 ^) is the Kolmogorov scale. For k < 1 this 



holds at least for K < z/" 1 £>! L" -1 . The energy balance equation ( |5.21| ), implies then 



that, for k > 1 , L^ 1 ^ K ^ rj' 1 , the mean energy flux 71-^ is almost equal to the mean 
injection rate ^><k and hence almost constant and equal to e in the inertial range. For 
k > 1 the same is true at least for L^ 1 <C K <C v~ x D\ L K ~ l . This confirms an intuitive 
picture of the non-dissipative energy cascade through the decreasing distance scales from 
the integral scale L where the energy is injected to the Kolmogorov scale rj at least for 
k > 1 . At the end the energy is dissipated on distances shorter than rj . What happens 
for k < 1 and K between 0[y~ x D\ L K_1 ) and 0(rj~ l ) remains to be understood. 



5.3 Gaussian case 



It will be instructive to look at the case when k = with 

D ij = 2D l 5 ij . (5.26) 

(Note that having finite D\ requires infinitesimal D when k — - > in order to renormalize 
the ultraviolet divergence in the momentum-space integral ( |4.3|) ; D will never show up 
below.) Although quite trivial, the case k = appears to be instructive and will serve 
below as the departure point for an expansion in powers of k. We immediately obtain 
for the k = operators : 

M S 2 C = 2D 1 V X1 ■ V X2 = -2D! A x , (5.27) 
MT = 2D 1 J2 V Xn -V x „, 

l<ra<n'<4 

= -2D! (A x + A Y + A z - V x • V Y - V Y ■ V z ) (5.28) 
in the difference variables X = xi — X2 , Y = x 2 — X3 , Z = X3 — X4 . The 2-point function 

F 2 (x) = (JUr^X*) = lkj e4k X l k r 2 ^(k) , (5.29) 

where 

C L (k) = J e~ ik ' x C(x/L) rf 3 x = L 3 C(Lk) . (5.30) 

Hence 

F a (x) = ^ L 2 /e ik ^ | k |-2c (k) . (5 . 31) 
The large L behavior may be obtained by expanding e* k ' X//i to the second order: 

^« = tk L 'f i k r 2 c » + ^,ikr 2 c» 

= ^ ((-AJ^OCO) + ^ (ftC-A^CXO) x l 

+ i757 ((Si^(-A)- 1 - ^-)C)(0)iV + 1^7^(0) |x| 2 + ^(L" 1 ) . (5.32) 



Thus, up to a constant and linear terms (the latter absent for real C) with C-dependent, 
L-diverging coefficients and a quadratic zero mode of Aif with a C-dependent coefficient, 
we obtain in the limit L — > oo the universal scaling solution F| c (x) = g^-|x| 2 , in 

accordance with the general result ( |5.9| ). Note that the separation of the scaling solution 
from the zero mode of the same homogeneity is ambiguous. 

The 4-point function in the stationary state is given by Eqs. ( |4.29| ) and (|4.30|) for 
M — 2, i.e. it may be written as a sum of contributions of three 2-particle channels: 

4 

< Y[ T(x n ) >= -F 4 (xi,X2,X3,x 4 ) + F 4 (xi,X3,X2,x 4 ) + F 4 (xi,x 4 ,x 2 ,x 3 ) (5.33) 

n=l 

with the single channel function 



F 4 = Ml 1 {M 2 X <g> 1 + 1 <g> M?) C L ®C L . (5.34) 
Working in the difference variables, we obtain for k = using the scaling operators (|5.27| ) . 



F 4 ( Xl ,x 2 ,x 3 ,x 4 ) = j^J e i ^+- +k ^ (£ kn-M -1 



n<rr 

((k! • k 2 )- x + (k 3 • k,))- 1 C L (kO C L (k 3 ) 
(27r) 3 <5(k 1 + k 2 ) (27r) 3 5(k 3 + k 4 ) f[0 

n=l 



1 I e i(k 1 .(x 1 -x 2 )+k3-(x 3 -X4) ( k 2 + k 2)-l ((fc^-l + ^2)-!^ 
4£> 1 J 



c L (kO c L (k 3 ) d3kl d3k3 



(2tt) 3 (2tt) 3 

= F 2 ( Xl -x 2 )F 2 (x 3 -x 4 ) , (5.35) 

see Eq. (|5.29| ). It follows from Eq. ( |5.33|) that for k = the equal time 4-point function is 
expressed by the 2-point function by the standard Gaussian formula. This remains true 
for the higher point functions. Let us consider the 6-point function 

6 

<]J T ( x n)>= Yl i7 6(x (J (i ) ,...,x (T(6) ) , (5.36) 

n=l aeS e /S 3 xS% 

where F§ is the contribution of a single 2-particle channel, 

^6 = M 6 l Yl ° <g> 1 2 ) (M2 1 ® 1 4 )) C L ® C L , (5.37) 

Explicitly, for n = , 

F 6 ( Xl ,...,X 6 ) = ^3 | e *(ki-(xi-x 2 )+k 3 .(x3-x 4 )+k 5 .(x 5 -x 6 ) ( k 2 + k 2 + k 2)-l 
■ ((k? + k^)- 1 + (k? + kg)" 1 + (kg + kg)" 1 ) ((k?)^ 

+(kir + (k^)- 1 ) &(ko c L (k 3 ) c L (k 5 ) 

= F 2 (x 4 - x 2 ) F 2 (x 3 - x 4 ) F 2 (x 5 - x 6 ) . (5.38) 
It should be clear how to extend this argument to a general 2M-point function. 



If we want to separate in F 2 m the non-universal and the universal terms, we have to 
extract from equation (|5.31|) more detailed knowledge about the 2-point function than 
that given by the expansion (|5.32|) . Since F 2 starts with an 0(L 2 ) term, we have to 
include into the expansion of, for example, _F 4 all terms of F2 down to the order 0(L~ 2 ) . 
Since 



*i(x) 



^Ei^a-AJ-ArcKoix"-,'" + o(l- 



n=0 



(5.39) 



F 4 (X,Y,Z) = ^ £ h. L 4 ~ m ~ n (d h . . .d im A _1 C)(0) 

m,n>0 
n + m<4 

■(d jl ...d jn A- 1 C)(0)X il ---X im Z h ---Z jn + 0{L- 1 ) . (5.40) 
For rotationally invariant C , this reduces to 

F 4 (X,Y,Z) = i l r (L 4 (A- 1 C)(0) 2 + iL 2 C(0)(A- 1 C)(0)(|X| 2 +|Z| 2 ) 

- m ( AC )(°) ( A ^)(0) (|X| 4 + |Z| 4 ) + i C(0) 2 |X| 2 |Z| 2 ) + G(L- 2 ) . (5.41) 

Note the presence of the universal Gaussian contribution F| C (X) F| C (Z) = |X| 2 |Z| 2 
on the right hand side. The rest of the terms are non-universal zero modes of the operator 

M S 2 C ® Mf {M S 2 C ® 1 + 1 ® Mf)~ x Mf (5.42) 

equal to Ai s 2 ®Ai s 2 = ADf AxAz in the action on Y-independent functions. That this 
must be the case may be easily seen by rewriting the relation ( |5.34j ) as the equation 

(M S 2 C ® M S 2 C ) (M S 2 C ® 1 + 1 ® M s 2 )~ l M S A C F 4 = C L ®C L . (5.43) 

Note that the separation of the universal Gaussian contribution F| C (X) F| C (Z) from the 
terms with m = n = 2 is ambiguous since there are non-universal terms of the same 
homogeneity degree. Nevertheless, the calculation for k = confirms the general struc- 
ture of the correlations with the non-universal terms given by homogeneous zero modes 
of differential operators. 



5.4 4-point and 6-point function for non-zero k 



Let us consider the 4-point function given by Eqs. ( |5.33|) and ([5.34 ) for general k between 
and 2. Recalling that F 2 = JA^Cl, it may be convenient to view Eq. ( 5.34 ) as the 
differential equation for F 4 : 



M A F 4 = F 2 ®C L + C L ®F 2 . 

The connected part F% = F 4 — F 2 <g> F 2 of F 4 satisfies the equation 

M 4 F< = C(F 2 ®F 2 ) , 

where C = (£) 1)3 + {K) lA + (£) 2(3 + (/C) 2j4 , see Eq. ( ggg ). 

First, let us study the L-dependence of F 2 ®F 2 . Since, as follows from Eq. 
v = 0, 



FoU) 



1 L 



2-K, 



\2-K 



3(2-k)Ui 



+ 0(L- 



(5.44) 

(5.45) 

) with 
(5.46) 



we infer that the product of 2-point functions 

F 2 (X)F 2 (Z) = 7 2 L 2 < 2 - K ) + 7 - I -^-L 2 - K (|X| 2 - K + |Z| 2 " K ) 

+ ^^Dl |XMZ| 2 -* + O(L-) . (5.47) 

Compare this expression with the similar expression ( [5 .41 ) for k = . The term ~ 
(|X| 4 + |Z| 4 ) present there comes from the 0(L~ K ) terms of the k > case. 

Let us return to the connected contribution to F 4 as given by Eq. ( |5.45|) . In the action 
on translationally invariant functions, 

Ma = - D^(X)d X id xi -p ij {Y)dYid Y i - D ij (Z) d Z i d zj 

- (D ij (X + Y) — D ij (X) - D ij (Y)) d x , d YJ 

- (p lj (Y + Z) - D lj (Y) - W{Z)) dyi d ZJ 
r D ij (X + Y + Z) - D ij (X + Y) - D ij (Y + Z) + D ij (Y)) d X i d Z] (5.48) 



for v = . Note that 

£(F 2 ®F 2 )(X,Y,Z) = ( D %3 (X + Y + Z) — D v (X + Y) 



-D^(Y + Z) + £>«(Y)J (^F 2 )(X) (djF 2 )(Z) (5.49) 

has a well defined limit when L — > oo since it involves only the derivatives of F 2 . In 
particular for m = and z/ = 0, the L — > oo limit of £ (F 2 ® F 2 ) is equal to 



9(2-re) 2 D 2 



and is a homogeneous (rotationally-invariant) function of X, Y, Z of degree 2 — k whose 
explicit form may be easily found. It is not difficult to write down a solution F£ sc of Eq. 
( |5.45|) for the limiting case. Indeed, for 

pcsc _ e 2 /|y|2(2-k) i |!7|2(2-k)\ £ 2 | Y |2-K|y|2->t fK KT\ 

^4 - 6(2 _ k) 2 (5 _ k)d 2 U A I +1^1 J 9(2-it) 2 r>? |A| |£| , (O.OIJ 

one obtains with the use of the decomposition A^ = .M 2 C (g) 1 + 1 ® .M 2 C — £ sc the 
required relation 

(^I c Ff)(X,Y,Z) = 9{2 f K?Dl C sc |X| 2 ^|Z| 2 - (5.52) 

since .M 2 c £*Dl + l<8>.M 2 c does not contribute to the left hand side and C sc annihilates 
functions depending only on X or only on Z . By the same argument as for the two-point 
function, the solution for finite but large L should differ from the universal scaling form 
by zero modes of Ai\ c so that 

Fl - E L 2 ^-^Y,lnmFl nm ^ Ft c (5.53) 

0<C4,n<2(2- K ) m 

where F£ nm are homogeneous zero modes of A4\ c of degree and the non-universal 
coefficients j nm depend on the source covariance C . In fact, a more complicated behavior 
with logarithmic terms 

Ft - E L^-^Y,ln mk (\n(Lf nm )) k Fl nm ^ Ft c (5.54) 

0<C4,«<2(2-k) m,k °° 



is possible if M.^ has zero modes (lnf nm ) 1 F nm , I < k, for f nm homogeneous functions 
of degree —1. 

The main problem which remains in the analysis of the universality of the 4-point 
functions is to find homogeneous zero modes of M. s £ with degree between zero and 2(2— k) 
and to decide which ones appear in the relation (|5.53|) . The study of such zero modes 
seems to lie beyond the scope of the techniques |[27|| developed for the singular elliptic 
operators of the type we encounter here. Note that in the full connected 4-point function 
only the zero modes with permutation symmetry survive. If the sources are not only 
homogeneous but also isotropic, only SO (3)-symmetric zero modes may appear. In fact, 
one should expect that such zero modes do enter into the expansions ( |5.53| ) or ( |5.54j ) 



for generic C whenever not forbidden by symmetries (which is a somewhat tautological 
statement). Existence of the zero-mode contributions to the 4-point function is important 
from the practical point of view. If we study a combination of 4-point functions (e. g. the 
4 th structure function) and we fix the range of variation of distances between the points, 
then the non-universal terms with the smallest CfT < 2(2 — k) dominate the behavior for 
large L and we would see the anomalous exponent C™™ rather than 2(2 — k) as governing 
the distance behavior. The possible appearance of zero modes of .M 4 C in the asymptotic 
behavior of the 4-point function would lead to the picture of restricted universality, 
holding strictly only for the 4-point correlator infrared-renormalized by subtracting 
homogeneous terms with non-universal coefficients typically divergent when the integral 
scale L grows to infinity. In the presence of non-universal terms with £ 4in = 2(2 — k) with 
are L- independent, the separation between the universal and non- universal terms in the 
4-point function would become ambiguous and the Eq. ( |5.51| ) would involve an arbitrary 
choice. 

Notice that adding the connected and disconnected contributions to _F 4 , one obtains 
(in the absence of logarithmic terms) 

F 4 (X, Y, Z) = G(X) + G(Z) + J2 L 2 ^-^ £ 7nm F 4 c inm (X, Y, Z) + o(L) (5.55) 

C4,n m 



for G(x) = i 7 2 £ 2(2 - K) + 7 3iri^£ 2 - K |x| 2 - K + 1H2 _J {5 _ K)D u |x| 2 ( 2 ->. Only F; 



4,mn 



terms may contribute to the structure function 5 4 (r) , as all the terms dependent on 
only one variable difference do not survive in S4. Hence the presence of the non-trivial 
contributions from the (generally) non-universal terms F£ nm is a necessary condition for 
the very non-triviality of S4 for large L . 

The 6 point function may be treated similarly. The connected part of the single 
2-particle channel contribution F 6 is given by 

F 6 (x a , . . • , x 6 ) = F 6 c (x!, . . . , x 6 ) + F 4 c (x 1; . . . , x 4 ) F 2 (x 5 , x 6 ) 

+ Fl (xi, x 2 , x 5 x 6 ) F 2 (x 3 , x 4 ) + Fl (x 3 , . . . , x 6 ) F 2 (x 1; x 2 ) 

+ F 2 (xi,x 2 ) F 2 (x 3 ,x 4 ) F 2 (x 5 ,x 6 ) , (5.56) 
see Eq. ( [4.3 1| ). It satisfies Eq. ( |4.34j ) which reads 



M 6 F< = ((£ 1 ,3 + £ 2 , 3 )(F 4 c (l,2)F 2 (3)) + (/: 1 , 2 + /: 2 ,3)((F 4 c (l,3)F 2 (2)) 

+ (£ 1 , 2 + £ 1 , 3 )(F 2 (1)F 4 C (2,3)), (5.57) 



where the numbers in the paranthesis label the 3 consecutive pairs of points. In the case 
that the non-universal non-constant zero modes do not appear in the connected 4-point 



function, the right hand side of Eq. ( p.57| ) has a finite L — > oo limit so that the relation 
of the type fl5.53|) , 



Fe - E L^~^Y.ln m Fl nm — > F 6 CSC , (5.58) 

0<Ce, n <3(2- K ) m 

with Fg mn being homogeneous zero modes of the degree £ 6jn of should hold for m = 
= v (with eventual additional logarithmic terms, as in (|5.54 ) ) . As we shall see in the 



next Section such zero modes with rotational and permutation symmetries exist at least 
for small k. Hence, the non-constant non-universal terms should appear if not already 
in the connected 4-point function then in the connected 6-point function. Whether they 
are visible in the corresponding structure functions is a more difficult question to which 
we do not have definite answer yet. 

Let us end this section by stressing a difference between the infrared renormaliza- 
tion required by the correlation functions of the PS and the renormalization of infrared 
divergences occurring in massless models of statistical mechanics or field theory. There 
the divergences are an artifact of the perturbative expansion and signal that the control 
parameters of the system are off their critical values. The divergences may be removed 
by finite shifts of the control parameters which depend non-smoothly on the coupling 
constants (hence the infinities in the perturbative treatment). On the contrary, in the PS 
model of turbulent advection, the genuine divergences appear when the size of the system 
tends to infinity and the correlators become finite only after the subtraction of the diverg- 
ing terms. Similar picture appeared in Polyakov's 2-dimensional conformal turbulence 
|p0|1 where it was argued that extracting the conformal-invariant part of the correlation 
functions requires subtraction of non-universal polynomial terms in the position spaceQ 
We expect the infrared renormalization of the type described above to be the general 
feature of universality in turbulent systems and the most important lesson to learn from 
the study of the PS even if, at present, we are unable to decide whether non-universal 
terms are present in the structure functions, the most commonly studied correlators. 

5.5 ^-expansion 

We may study the homogeneous zero modes of the operator M. s £ in perturbation expan- 
sion in powers of k. We shall work here only to the first order in k. Eq. fl4.6p implies 
that (for m — ) 

D ij (x) = 2D 1 5 ij + kD x (25 ij In |x| - xV|x|~ 2 + 5 ij ) + 0(k 2 ) 

= 2D 1 5 ij + 2KD 1 R ij (x) + 0(k 2 ) . (5.59) 

As a result, to the first order in k, 

M S 4 C = Mf >0 - 2KD 1 R i i{X)d X id xi - 2kD x R l] (Y) d Y , d Yi - 2kD x R ij (Z) d Z i d zj 

- 2kD 1 (i? lJ (X + Y) — R ij (X) - R t3 {Y)) d X i d Y3 

- 2kD 1 (r 11 (Y + Z) - R lJ (Y) - R ij (Z)) d Y r d Z3 

- 2kD 1 (i#'(X + Y + Z) - R ij (X + Y) - R ij (Y + Z) + R ij (Y)) d xi d ZJ 

= Mf )0 + 2k D l Vi , (5.60) 



3 K.G. thanks A. Polyakov for discussion of this point 



where the subscript "0" refers to k = 0. M% commutes with 3- dimensional translations 
and rotations and with the permutations of four points. We shall search for its zero modes 
respecting these symmetries. Note that 

M% Q = -2D X (A x + Ay + Ag) , (5.61) 

where 

X = X, Y = a/2(Y + ~X + ~Z) , Z = Z. (5.62) 



Denoting i? = (X 2 + Y 2 + Z 2 ) 1 / 2 , we obtain 

Ml, = - 2 -§^d R R 8 d R + ^$ (5.63) 



where <3> is the Laplacian on the sphere S 8 in the space of (X, Y, Z) . The symmetric 
zero modes of the lowest homogeneity of M.^ Q occur in degree zero (constants) and in 
degree 4 . The latter, when expressed in terms of point-differences, have the form 

a^(x n -x n /) 4 + 6^(x n -x m ) 2 (x n -x m /) 2 + c^(x n -x n /) 2 (x m -x m /) 2 , (5.64) 

{n,n'} {{n,m},{n,m'}} {{n,n'},{m,m'J} 

where the pairs {n,n'} and {m,m'} are assumed different, as well as the pairs {n,m} 
and {n,m'} and where 

20a + 286 + 6c = . (5.65) 

The constant survives as the eigenvalue of A^| c for k ^ 0. The change of the two 
independent zero modes R A fi , i — 1, 2 , of the form (|5.64|) may be traced by the standard 
degenerate perturbation theory. This is done as follows. Looking for a homogeneous zero 
mode R 4+kX {a\ fx + a 2 f 2 + k / 3 ) with a homogeneous degree zero function f 3 orthogonal 
to fi t 2 in L 2 (S 8 ) , we obtain in the linear order in k 

A^^ (A J R 4 lni?(a 1 / 1 + a 2 / 2 ) + R 4 f 3 ) + 2D 1 V,R 4 (aj, + a 2 f 2 ) = (5.66) 

or, using the form ( |5.63| ) of Ai s ^ , 

-15A(ai/i + o 2 / 3 ) - 44/ 3 + $/ 3 + 4^4^(01/1 + 02/2) = . (5.67) 



Upon taking the L (S ) scalar products with /i^, f 3 drops out resulting in the relation 
-15A ^ (/,../• • ^V 4 R^fj) aj = (5.68) 

j=l,2 i=l,2 

Hence A has to solve the equation 

det ^Vt&fj) - A(/i, /,•)) = . (5.69) 

It may be decided numerically whether this equation has solutions A < —2 so that there 
exist homogeneous zero mode of M. s £ with degree < 4 — 2k for small k which, multiplied 
by 7 L-( 2+A ) K , should appear as non-universal terms in the truncated 4-point function. 
Given A and a non-zero solution (ai,a 2 ) of Eq. ( |5.68[) , the O(k) contribution / 3 to the 
zero mode should be computed from Eq. ( |5.67| ). 



If Eq. ( |5.69p ha s one double root but there is only one solution (01,02) (up to nor- 
malization) of Eq. fl5.68|) then one should look for an additional zero mode of A4| c of the 
form 

k (InR) R 4+KX { ai h + a 2 f 2 + k/ 3 ) + R A+KX {hh + b 2 f 2 + «/ 4 ) . (5.70) 
This leads in the first order to the relation 

- 15 ((ai + Xh)fi + (aa + A6a)/ a ) - 44 f 4 + $/ 4 + ^^fa/i + b 2 f 2 ) = (5.71) 
which implies that 

- 15 ]T (fi, f 3 ) ( aj + Aft,-) + ]T (/ f , ^ V.R'fj) b 3 = . (5.72) 

j=l,2 j=l,2 

The latter equation has a non-zero solution for (6 1; 6 2 ) • In this case, the logarithmic terms 
~ (InL) L _ ( 2+A ^ K should appear in the asymptotic expansion of F£ for k > if A < —2. 

The zero modes of the scaling operator Ai S Q may be studied similarly by writing 

Mf = M% +2kD 1 V 6 , (5.73) 

and using the perturbation expansion. By a linear change of variables, -Mg c is equivalent 
to the Laplacian in R 15 . As for -M| c , the homogeneous zero modes of A^| c with 
the lowest degree, symmetric under rotations and permutations, are constants and the 
4 th -order polynomials ( |5.64j ) satisfying this time the relation 

20a + 566 + 36c = . (5.74) 

For k > the constants are still annihilated by .Mg c whereas the two 4 th -order zero 
modes may change slightly the homogeneity degree (one of them may also pick a log- 
arithm). As for small k the deformed degree will still be smaller than 3(2 — k) , the 
deformed zero modes, accompanied by a positive power of L and a non-universal co- 
efficient, should appear in the asymptotic expansion of the connected 6-point function 
discussed above. 



6 Conclusions 



At the beginning of these lectures, we have sketched the main idea of Kolmogorov's pic- 
ture of universality in the fully developed turbulence, related to the energy cascade. We 
briefly described how the turbulence problem may be formulated in the language formally 
very close to that of equilibrium statistical mechanics or Euclidean field theory stressing, 
however, the important differences between the two classes of problems. The bulk of the 
lectures was devoted to the discussion of a simple model of turbulent phenomenon, the 
passive scalar (PS), which we treated as a test ground for the universality ideas. The 
exact solution for the model in terms of singular elliptic many-body operators showed 
non-universal diffusive decay of correlations for deterministic initial data but a universal 
super- diffusive decay for random homogeneous initial data. In the presence of random 
sources injecting the energy at small wavenumbers, the steady state exhibited an inertial 
range with the energy cascade and universal scaling in the connected equal-time two- 
point function of the scalar. We have argued that the higher-point equal-time connected 



functions are expected to show only restricted universality with the terms dominating 
for large size systems characterized by universal anomalous exponents and non-universal 
amplitudes. We have shown that such terms should appear in the 6-point function if not 
in the 4-point one, at least for small velocity exponent k . Further work along well traced 
direction is required to decide on this question as well as on the question what terms dom- 
inate the respective structure functions. One of the open problems is to relate the terms 
with anomalous exponents to the physical idea of intermittency. One should also under- 
stand what role, if any, is played in the model by the volume preserving diffeomorphisms. 
Upon infrared renormalization which subtracts the non-universal terms, the equal-time 
correlators of the scalar should become finite and universal (with normal scaling) in the 
limit when molecular diffusivity v — > and the integral scale L — > oo. In particular, 
no //-dependent subtractions are needed in the steady-state equal-time correlators of the 
scalar and their short-distance behavior is easy to analyze JI4] . This behavior de- 



termines what ultraviolet (short distance) renormalization is required by the composite 
operators involving the gradients of T . The basic role of v is to assure the convergence 
to the steady state by providing the mechanism which removes on small scales the energy 
injected by the sources at large distances. 

We expect the main features of the behavior of the PS described above to be present 
also in other systems exhibiting a fully developed turbulence, including the Navier-Stokes 
flow. Still much remains to be done even for the PS. One important step could be a 
complete perturbative analysis with the above picture of restricted universality established 
order by order. The k expansion discussed above is a possibility but already the analysis 
of the lowest terms proved to be quite complicated p6[ . Other possibilities are the quasi- 
Lagrangian perturbation expansion or the expansion proposed in ||14| . Ultimately, if no 
way to compute exactly the non-universal terms in the correlators is found, a multiscale 
analysis employing perturbative arguments in conjunction with a renormalization group 
type analysis should provide a tool to fully control the behavior of correlation functions. 
Developing such methods for the PS may be a first step towards a complete theory of fully 
developed turbulence. More than half-century after Kolmogorov's work, despite further 
progress, a fully developed understanding of the fully developed turbulence remains a 
major challenge. 
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